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Users  guide  to  SROSOL 
Introduction 

This  report  briefly  describes  a  computer  code  called  SROSOL  written  by  one  of  the 
authors,  Paul  S.  Julienne,  that  calculates  wave  functions  and  matrix  elements  for  diatomic 
molecules.  These  matrix  elements,  for  example  Franck-Condon  factors,  can  then  be  used 
to  investigate  phenonema  such  as  molecular  radiative  emission  and  absorption. 

The  focus  of  this  report  is  to  serve  as  a  users  manual  for  SROSOL;  hence  the  primary 
emphasis  is  on  the  input  to  the  code  and  how  it  is  used.  The  input  parameters  are  discussed 
in  some  detail,  especially  for  a  Morse  potential.  Then  comparisons  of  results  from  this  code 
for  Franck-Condon  factors  for  the  First  Positive  Band  of  N2  are  given  with  two  similar 
calculations. 

Finally,  an  example  of  the  use  of  SROSOL  results  is  presented  in  the  form  of  a  calcu¬ 
lated  molecular  emission  spectrum. 

The  program  SROSOL  calculates  eigenvalues  and  normalized  bound  or  continuum 
solutions  to  the  one-dimensional  Schrodinger  equation 

d2i/>/dR2  +  (Att2 /h2)2fj,(E  -  U(R))il>  =  0 

where  /x  is  the  reduced  mass,  and  h  is  Planck’s  constant.  These  wave  functions  can  be 
used  to  calculate  matrix  elements  of  the  type 


<  ^u\0n\^2k  > 

where  is  the  wave  function  of  level  i  in  potential  1,  and  ip2k  is  the  wave  function  of  level 
k  in  potential  2:  6 12  is  the  operator  being  evaluated.  The  numerical  algorithm  is  adapted 
from  Gordon1. 

The  first  wave  function  is  calculated  on  a  self-determined  point  grid  which  optimizes 
the  propagation  intervals  both  for  the  self-propagation  and,  by  virtue  of  a  comparison 
with  the  second  potential,  for  the  calculation  of  the  matrix  element  later.  The  second 
wave  function  is  calculated  on  the  fixed  point  grid  established  by  the  propagation  of  the 
first  wave  function. 

The  actual  equation  solved  by  the  code  is 


d2j>/dx2  +  ((E  -  lf(R))/Eo)i>  =  0 


where 


E0  =  {h2/ Air2)/ 2^1 
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gives  the  working  energy  unit,  and  x0  gives  the  working  length  unit,  e.g.,  if  E0  is  in  cm  1 
,  x0  is  in  Bohr(Bohr  radii),  and  fi  is  in  atomic  mass  units(amu),  then 

E0  =  60.19974//t 

in  units  of  cm-1.  The  energies  and  potentials  E  and  U  must  be  given  in  the  same  units 
as  Eo  and  x0.  The  choices  available  to  the  user  are  determined  by  the  parameters  IE,  IX, 
and  IM  defined  on  Line  3  of  the  input  file;  these  are  discussed  in  the  next  section. 

Although  SROSOL  can  treat  potentials  U(R)  which  are  Morse  type,  exponential  re¬ 
pulsive  type,  and  numerical  splines,  this  report  is  limited  primarily  to  Morse  potentials. 

Input  File 

This  section  is  a  line-by-line  discussion  and  description  of  the  input  file  FT05. 

The  line  numbers  appear  in  italics  followed  by  the  input  parameters  on  that  line.  A 
description  of  these  parameters  follows  each  line  number. 

Line  1 :  PRBLBL  (This  is  a  label  describing  the  problem  being  treated.) 

Line  2:  IPROB,  NPOT  (see  table  below) 

(value  of  IPROB:  wave  functions  conditions) 

0:  calculate  eigenvalues  only  (no  matrix  elements) 

1:  ipu  and  ip2k  3X6  both  bound  states 

2:  ipu  is  bound  and  ip2\ b  is  continuum  (any  energy) 

3:  ip\ i  is  bound  and  ip2k  is  a  continuum  state  of  the  same  energy  (i.e.,  predissocation) 
4:  ipu  and  ip 2k  are  both  continuum  states 
5:  unused 

6:  only  tabulate  and  graph  potentials  -  no  eigenvalues  or  matrix  elements 
7:  ipu  and  ip2k  are  both  continuum  states;  all  i,k  pairs  calculated 
8:  ipu  and  ip2k  are  both  bound  states;  only  i=k  pairs  calculated 


NPOT=  the  number  of  potentials.  The  range  is  from  2  to  5  when  IPROB  is  greater 
than  zero. 

Line  3:  IE,IX,IM,/z  (see  table  below) 


value 

IE 

IX 

IM 

1 

au(Hartrees) 

Bohr 

electron  masses 

2 

1/cm 

Angstrom 

amu 

3 

4 

eV 

reduced 

reduced 

reduced 

2 


fi  is  the  reduced  mass  of  the  molecule.  For  reduced  units,  Eo  =  1/2/x 

Line  4-  DELH,HO,HMIN,DELNOR  (see  table  below) 

DELH:  basic  numerical  accuracy  (0.001  generally). 

HO:  initial  propagation  step  (normally  set  equal  to  DELH). 

HMIN:  minimum  propagation  step  size  (normally  set  equal  to  DELH). 

DELNOR:  convergence  criterion  for  phase  of  continuum  wave  function  (generally  .001 

Lines  5-8  define  the  first  potential,  and  are  illustrated  for  the  case  of  a  Morse  potential. 
The  potential  is  read  by  subroutine  RDPOT,  and  evaluated  by  function  UII. 

Line  5:  LABEL1  (Label  for  the  first  state.  ) 

Line  6:  4,JROT,  where  JROT  =  rotational  quantum  number. 

Line  7:  1,ITYPE,4,  (ITYPE  is  described  in  the  following  table.) 

Line  8:  A1,A2,A3,A4  (see  table  below.) 


ITYPE 

A1 

A2 

A3 

A4 

1 

wc 

Dc 

Re 

Te 

2 

U)eXe 

Re 

Te 

3 

wexe 

De 

Re 

Te 

4 

fie 

Dc 

Re 

Te 

For  the  Morse2  potential, 

U(R  -  Re)  =  De(  1  -  e-KR-R^Y  +  Te, 

where 

De  =u>l/wexe, 

and 

/?c  =  1.2177X10Va/^Z^T- 

Here  De,  we,  a-nd  fie  axe  in  cm-1. 

For  the  exponential  repulsive  potential,  Line  6  becomes  ’5, JROT’,  and  ITYPE  =0  on 
Line  7.  The  form  of  the  potential  is  then 

U(x)  =  A2e~A3('x~Al)/A7  +  AA 

Lines  9  through  12  describe  the  second  state  and  are  similar  to  Lines  5  through  8. 

Line  13:  ITAB(l-5),IPLOT(l-5),XI,XF,XINC,EMIN,EMAX 

ITAB(I):  (=0,  do  not  tabulate  potential  I;  =1,  tabulate  potential  I.) 
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IPLOT(I):  (=0,  do  not  plot  potential  I;  =1,  plot  potential  I.) 

Tabulation  and  plotting  begin  at  R  =  XI  and  run  to  R  =  XF  in  increments  of  XINC. 
EMIN  and  EMAX  are  the  mininum  and  maximum  energies  for  the  plots. 

Line  14-  ITYPE 
Line  15:  FAC(I) 


ITYPE  Function 

0  FAC 

1  B*FAC 

2  FAC(1)+FAC(2) 

3  FAC(1)*EXP(-FAC(2)*R) 

4  Spline  f*FAC 

5  SpHne  f*FAC*B 

6  Spline  f*FAC2 


Constant 

Eq/R2,  rotational  function 
Linear 
Exponential 
Spline  function 

Spline  function  X  rotational  function 
Spline  function  squared 


For  example,  for  ITYPE=2,  FAC(1)=1.0,  FAC(2)=0  calculates  Franck-Condon  fac¬ 
tors. 

For  ITYPE=2,  FAC(1)=0,  FAC(2)=1.0  calculates  matrix  elements  of  R. 

Line  16:  IREAD 

IREAD=1:  read  eigenvalue  hst  for  bound  eigenvalues. 

IREAD=2:  read  energy  list  for  continuum  energies. 

IREAD =3:  search  for  eigenvalues. 

Line  17:  TXL,X0,TXR 

X0  is  approximately  equal  to  Re  for  bound  wave  functions.  TXL  and  TXR 
are  left  and  right  guesses  for  classical  turning  points.  For  continuum 
functions  when  IPROB=4  or  7,  X0  is  the  distance  up  to  which  the  wave- 
function  is  stored  (outer  range  of  integration).  When  IPROB=2  or  3,  X0 
is  the  distance  beyond  which  convergence  of  the  phase  is  checked  (also 
true  for  IPROB=4  or  7). 

Line  18:  NE,EI,EF,DELE 

NE  =  number  of  equal  intervals  between  El  and  EF  to  be  used  in  the  eigen¬ 
value  search,  the  interval  should  be  less  than  the  smallest  difference 
between  adjacent  eigenvalues.  Subroutine  EIGEN  determines  if  there  is 
an  eigenvalue  in  each  interval  and,  if  so,  converges  to  it. 

El  =  lowest  energy  in  eigenvalue  search.  It  should  be  larger  than  the 
potential  minimum  to  prevent  convergence  failure;  the  point  (X0,EI) 
should  be  in  the  classical  region  of  the  potential.  Beware  when  chang¬ 
ing  rotational  quantum  number  J. 


4 


EF  =  highest  energy  in  eigenvalue  search.  It  should  be  less  than  the 
asymptotic  energy  to  prevent  convergence  failure. 

DELE  =  convergence  parameter  in  eigenvalue  search.  The  eigenvalue  is  located 
to  within  DELE  of  its  actual  value  (more  correctly,  the  zero  of  the 
mixed  Wronskian  of  the  left  and  right  numerically  propagated  wave- 
functions  is  located  to  within  DELE  of  the  actual  zero  for  the  given 
propagation  parameter).  An  extra  linear  interpolation  usually  locates 
the  eigenvalue  to  within  approximately  0.1DELE.  Thus,  DELE  =  1  cm  1 
usually  results  in  eigenvalues  to  a  few  hundredths  of  one  cm  1. 

Lines  19  through  21  describe  the  second  state  and  are  similar  to  Lines  16  through  18. 

Applications 

Examples  of  the  use  of  SROSOL  are  given  in  the  following  pages.  In  Table  1,  a  sample 
input  file  is  shown  to  illustrate  the  descriptions  in  the  preceding  section.  In  Table  2,  a 
calculation  of  Franck-Condon  factors  for  the  First  Positive  band  of  N2  is  displayed,  along 
with  two  similar  results  from  the  literature;  agreement  is  good. 

An  estimated  emission  spectrum  for  this  same  N2  band  is  plotted  in  Figure  1.  For 
this  calculation,  the  electronic,  vibrational,  and  rotational  temperatures  were  set  at  5000° 
Kelvin.  Populations  for  the  various  vibrational  and  rotational  levels  were  obtained  from 
an  LTE  calculation5  .  A  pressure  of  100  atmospheres  was  assumed.  A  Voigt  profile 
was  assumed  for  each  fine  with  collisional  broadening  of  10  MHz/Torr  of  gas  pressure  as 
suggested  by  Steinfeld6  (see  also  Wang  and  Xia7 ). 

Figure  2  shows  an  emission  spectrum  of  the  First  Positive  band  of  N2  obtained  from 
a  discharge  tube  as  given  by  Keck,  et.  al8.  Virtually  no  information  about  this  spectrum 
is  given  in  Ref.  8;  however,  it  was  used  as  a  model  for  generating  the  calculated  spectrum 
shown  in  Figure  1.  Parameters  were  adjusted  in  this  calculation  to  try  to  obtain  essentially 
the  same  spectral  shape  given  in  Figure  2.  It  is  not  known  if  the  assumed  pressure  and 
temperatures  are  correct,  or  if  LTE  applies  in  this  case. 

For  a  comprehensive  treatment  of  N2  molecular  emission  and  absorption,  several  bands 
must  be  included  in  the  spectra.  It  would  be  necessary  to  apply  SROSOL  to  each  band 
to  obtain  Franck-Condon  factors  and  matrix  elements  in  order  to  be  able  to  calculate  the 
contribution  of  each  band.  Such  calculations  have  been  attempted  for  N2  and  S2. 
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Table  1:  Sample  SROSOL  input  file 


Line  1  :  N2  1st  Positive 

Line  2:12 

Line  3:2  2  2  7.00335 

Line  4  :  0.001  0.001  0.001  0.0000001 

Line  5  :  N2  B  3  Pi  g  State 

Line  6:40 

Line  7:124 

Line  8  :  1733.39  14.122  1.212600  59619.4 
Line  9  :  N2  A  3  Sigma+u  State 
Line  10:  4  0 
Line  11:  12  4 

Line  12:  1460.64  13.872  1.286600  50203.6 
Line  13:  1100000000  0.5  3  .1  0  10000 
Line  14:  2 
Line  15:  1.0  0.0 
Line  16:  3 

Line  17:  0.8  1.212  2.0 

Line  18:  5  60000  90000  0.000001 

Line  19:  3 

Line  20:  0.8  1.293  2.0 

Line  21:  5  50300  75000  0.000001 

This  input  file  is  constructed  to  calculate  Franck-Condon  factors  for  the  First  Positive 
band  of  Nitrogen.  On  Line  14,  ITYPE  =  2;  on  Line  15,  FAC(l)  =  1.0,  and  FAC(2)  =  0.0. 
NPOT  =  2  (Line  2).  On  Line  3,  IE  =  DC  =  IM  =  2,  which  gives  energy  in  l/cm,  length 
in  Angstroms,  and  mass  in  amu.  On  Lines  7  and  11,  ITYPE  =  2;  on  Lines  8  and  12,  Al, 
A2,  A3,  and  A4  are  defined  by  the  following  equations: 


Al  =  we, 

A2  =  ojexe, 

A3  =  Re, 

A4  =  Te. 

On  Line  16,  IREAD  =  3,  which  tells  SROSOL  to  search  for  eigenvalues.  The  param¬ 
eters  on  Lines  17  and  18,  and  Lines  20  and  21,  are  obtained  from  potential  plots  U(R)  of 
the  states  involved  in  the  calculation. 
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Table  2;  comparison  of  Franck-condon  factors  from  3R050L  with  other  calculations 
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N2  (1st  Positive)*  B  3  Pi  g->A  3  Si+u 
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Figure  1.  Calculated  emission  spectrum  from  the  First  Positive  Band  of  N2. 
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Figure  2.  Discharge  tube  spectrum  for  N2  from  Keck,  et.  al8. 
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